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Abstract: We construct a Lax operator for the G2-Calogero-Moser model by means of 
a double reduction procedure. In the first reduction step we reduce the Ag-model to a 
i?3-model with the help of an embedding of the -£?3-root system into the Ag-root system 
together with the specification of certain coupling constants. The G2-Lax operator is 
obtained thereafter by means of an additional reduction by exploiting the embedding of 
the G2-system into the -E^-system. The degree of algebraically independent and non- 
vanishing charges is found to be equal to the degrees of the corresponding Lie algebra. 



1. Introduction 



^ ■ The Calogero-Moser models ||, |||, || ^, 0, § constitute a large class of well studied 



interacting many particle systems. The models are very universal in the sense that they 
can be cast into a form in which the potential term includes a sum over all roots a of some 
root system A and the functional dependence of the potential is V(x) ~ 1/ sn 2 (x), with sn 
being an elliptic function together with its various limits l/sinh 2 (x), 1/ sin 2 (x) and 1/x 2 . 
Often it is useful to treat the latter cases independently for their own sake. Due to their 
universal nature the models find a wide range of applications in physics, as for instance 
to characterize anyons on the lowest Landau level ||, to describe certain properties of 
quantum Hall droplets |]l0| and in various ways in conformal 12, 13, [l4j] and boundary 



1 15] conformal field theories. 

The Hamiltonian for an n-particle Calogero-Moser system reads 



n = ^-\Y,aly^-<l) 9a^iR, q, P eR n , (1.1) 



with n being the dimensionality of the space in which the roots a are realized. At this point 
we impose only the restrictions g a = g- a on the coupling constant, even though later on we 
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equate more of them for reasons to be explained. One of the most prominent feature of these 
models is their integrability, meaning here the existence of a sufficient number of conserved 
quantities (integrals of motion) Ik in involution. A standard technique to construct these 
charges, the so-called isospectral deformation method, goes back almost forty years [16]. It 



consists of formulating Lax pair operators L and M as functions of the dynamical variables 
qi and pi for 1 < i < n, which satisfy the Lax equation L = [L,M], upon the validity of 
the classical equation of motion resulting from Ql.l| ). The Lax operator is then the starting 
point for the construction of conserved charges of the form 1^ = tr(L k )/k, with I2 ~ TL, of 



classical r-matrices [JITJ, |18|, |19|, |20|, |2l|, |22J], spectral curves |23|, |24], |25|, |26|] and various other 
important quantities. For root systems A, which can be associated with a Lie algebra g, 
i.e. crystallographic ones 1 , a natural Ansatz is to expand L and M in terms of the elements 
H,E a of g 

L=p-H+^2g a f{a-q)E a and M = m ■ H +^ g a h(a ■ q)E a . (1.2) 

Alternatively, one may also expand L and M in terms of other non-commuting objects, 
such as Coxeter transformations, and perform a similar analysis p8| , 29]. Substitution of 
these operators into the Lax equation yields various constraining equations, which for a 
given potential determine the functions f(x) and h(x) in L and M as defined in equation 
(|l.2| ). We choose here as convention the Cartan-Weyl basis tr(HiHj) = 5ij, tr(E a E^ a ) = 1, 
which is consistent with the well-known commutation relations (e.g. |3C| ]) 

[Hi,Hj] = 0, [Hi,E a ] = a l E a , [E a ,E- a ] = a - H, [E a ,Ep] = £ aj pE a+ p. (1.3) 

Then by direct substitution it follows that the Lax equation holds once the functions 
fa{x) = g a f(x) and h Q (x) = g a h(x) satisfy 

g(x) = f'(x), p = ^2af a (a-q)h- a (-a-q), j-m= ^ e«,/3 ^— 7 ^-^y q \ (1.4) 

0+13=7 



Assuming further that I2 = T~C and the classical equation of motions resulting from (1.1), 
one obtains three additional equations 

f( x )f(- x ) = -V(x), p=-^J2»9 2 aV'(a-q) and f(x) = ±f(-x). (1.5) 

For the stated potentials it is straightforward to use the factorizing condition in ( |1.5| ) to 
determine f(x), i.e. 1/ sn(x), 1/ sinh(x), 1/ sin(x) and 1/x, and therefore g(x) by taking its 
derivative. Thus to establish the integrability of the system ( |l.l|) reduces to the question of 
whether the third set of equations in (|1.4j ) admits a solution for the vector m and therefore 
guarantees the existence of the operator M. As these equations are in general highly 



1 For non-crystallographic root systems one may exploit the fact that they are embedded into crystallo- 
graphic ones and use a reduction procedure to obtain a meaningful Lie algebraic operator J27|. 
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overdetermined, the answer to this question depends crucially on the structure of the Lie 
algebra. As was observed long time ago |31], for given potentials as in ( |1.1[ ) the relation 



f(x)f'{y) - f'(x)f(y) = f(x + y) [V(x) - V(y)} 



(1-6) 



holds, such that the equations may be simplified further. Then the last set of equations in 
fll.4| ) reduces to 

57(7 ■m)= ^2 e a,p9a9/3V(a-q) = ^2e an g a g a+y V(a-q). (1.7) 

a+/3=7 

Clearly 

"^E^,^'^, (1.8) 

a6A i=l 

with Aj being a fundamental weight, is a solution to ( |1.7|) when 7, is taken to be a simple 
root. However, for (|ll^) to be a proper solution one also has to verify whether it solves 
Ql.7| ) for the remaining roots. In summary, we can say that if the system ( |1.7[ ) can be 
solved for the vector m for a particular Lie algebra g, then the system ( |1 . 1| ) is classically 
integrable. The reverse statement does not hold. 



2. The G 2 -Lax operator 

It turns out that only when the algebra g in (|1.2| ) is taken to be Ag one obtains directly, 
meaning that all quantities in ( |1.1| ) belong to Ag, a solution for the Lax operator with 
the condition that the corresponding equation of motion holds. In all other cases one 
needs to device alternative methods. The expressions for the Bi, Ci and .D^-algebras were 
obtained |32], |?|, || from ^^-theories by specific transformations of the dynamical variables 
and a subsequent constraint on certain coupling constants g a . For the remaining algebras 
different types of techniques have been developed [|32], [7|, ||, ^6], ^8], |29| . Surprisingly for the 
Lie algebras -£/6,7,8> -^4 an d G2 no Lax pair was known until fairly recent [26]. 

In particular, the latter, the G2-Calogero-Moser model, constitutes a standard simple 
example, since it can be viewed as the classical three-body problem with a two and a three- 
body interaction term [ 33 1 . For a specific realization of the roots (see below), the potential 
term in the rational case simply reads 

i<j<i<3 VH H3) i<«<i<3 yH H ° H ' 

with g s , gi being coupling constants. It appears to be rather surprising that despite the 
simplicity of this model, apart from the expressions in |26]| , a general simple formula for 
the Lax operator along the line of the original work of [^2] , 0, ||] may not be found in the 
literature. It will be the purpose of this paper to provide such a simple expression. 
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2.1 Direct computation 

Let us commence by directly analyzing equation ( |1.7| ) for G2. For the explicit calculation 
we require first the roots of G2. We recall the general fact, see e.g. [|34], [35], |36|, that the 
entire root system can be generated by h — 1 successive actions of the Coxeter element a 
on bi-coloured simple roots, i.e. 7 = ± cti, with h being the Coxeter number. It turns 
out to be convenient to abbreviate the roots accordingly, that is we define a p 7j 
1 < i < £ =rank g and < p < h — 1. 

For G2 we have I = 2, h = 6 and the 2 x 6 = 12 roots are computed to 



a itP for 



i\p 





1 


2 


3 


4 


5 


1 


— OL\ 


-(2&1 + a 2 ) 


— (&i + a 2 ) 




2ai + a 2 


OL\ + Q!2 


2 


OL2 


— (3ai + 02) 


-(3ai + 2a 2 ) 


— OL2 


3ai + 02 


3ai + 2a2 



The roots a.i :P of the AG 2 -root system. 



In addition, we require the structure constants e a ^ for the analysis of (|T 

.2 

'a,/3 



The square 

of the latter can be fixed by means of the well-known formula e 2 g = a 2 n(m + l)/2, where 



the integers n, m are determined by the so-called a-string through /?, i.e. the largest values 

ma are still roots, see e.g. [37]. The overall signs are 



for n, m such that (3 + na and (3 
in general not fixed and are subject to convention. However, some consistency relations 
have to hold, resulting from the anti-symmetry of the commutator, the reality condition 
and the Jacobi identity when a + (3 = 7 



£a,f3 — ^/3,« — ^/3,— 7 — 7,/3 — E-a,— 7 — £—7,0 — (^-2) 
£ — o,— (3 = £—f3,—a = /3,7 = £7,— /3 = £—0,7 = £7,-0- (2-3) 

We choose here the short roots to have length a 2 = 2 and the long roots a2,p to have 
length a 2 = 6. As complete lists of structure constants are difficult to find in the literature, 
we present here a consistent choice for the 12 x 12 = 144 structure constants, with 60 of 
them non-vanishing 



"i,p\" i j',q 


"1,0 


"1,1 


"1,2 


"1,3 


"1,4 


"1,5 


"2,0 


"2,1 


"2,2 


"2,3 


"2,4 


"2,5 


"1,0 





M 


-2 





2 


A* 













-A* 





"1,1 


-M 





M 


2 





-2 














M 


-A* 


"1,2 


2 


-A* 





M 


-2 





-A* 
















"1,3 





-2 


-A* 







2 


A* 


A* 














"1,4 


-2 





2 


A 4 





-A» 





-A» 


A* 











"1,5 


-A* 


2 





-2 













-/^ 


A* 








"2,0 








M 


-A* 














-M 





A* 





"2,1 











-M 


A* 














A* 





-A* 


"2,2 














-M 


1* 


M 











A* 





"2,3 





















-A* 











A* 


"2,4 




-A* 














-A* 





A* 











"2,5 





A* 
















A* 





-A* 









The G^-structure constants £i, P ,j,q with = y/3. 
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In order to obtain the previous table we only fixed five signs by convention and deter- 
mined the remaining ones by means of the relations ( |2,2| ) and ($2^) . 

Having assembled all necessary data, we can present a simple argument which demon- 



strates that it is not possible to solve the constraint (1.7) directly. Choosing all coupling 
constants g^ to be either g s or gi for a to be a short or long root, respectively, we may for 
instance add up the equation ( |L7] ) for the three choices 7 = Si,o, 7 = ai,2 and 7 = 5i,4 
and find after cancellation of g s 

(51,0 + 5i, 2 + 5i )4 ) ■ m = = 2\/3gi [V(5i )0 ■ q) - V(a 2fi ■ q) + V(a 1A ■ q) - V{a 2A ■ q)] ■ 

Clearly, the right hand side is not zero in general, and hence we can not solve the constraint 
Ql.7| ) directly. One reaches the same conclusion by taking the expression in ( |1.&| ) and trying 



to verify (1.7) for 7 to be a non-simple root. 

However, if we switch off the two particle interaction, i.e. we take g s = in (|2.1| ), we 
may construct a particular solution. Taking for instance 5i = £\ — £ 2 , ol 2 = — 2ei + £2+^2 
as concrete realization for the simple roots of G 2 in M 3 , with £$•£,• = Sij, and setting g s = 0, 



we can solve (1.7) by 



m 1 = 0, (2.4) 

m 2 = ^= [V(a 2 , ■ q) + V(a 2A ■ q) - 2V(a 2 , 5 ■ q)} , (2.5) 

m 3 = [2V{a 2A ■ q) - V(a 2fi ■ q) - V(a 2 , 5 ■ q)\ . (2.6) 

To find a general Lax operator which involves all terms of the potential one needs to device 
other techniques. 

2.2 The G2-Lax operator from double reduction 

The construction procedure is summarized by the following Dynkin diagrams: 



ac\ a 2 03 04 05 «6 ^ -j- 6l\ a 2 \ 03 



In the first step we start by folding the Ag -root system to the A^^-xoot system by means 
of a map oj. Subsequently we constrain some of the coupling constants through a map r, 
which amounts to an elimination of some particular roots of Aq, such that we obtain two 
copies of a i?3-root system. From this system we obtain the G2-root system, by the action 
of a further map a). 

The L-operator for the G2-Calogero-Moser model then reads 

6 6 

L(p, q) = uj^uj^p ■ H + T (9i, P )f{uujai :P ■ q)E aip (2.7) 

i=l p=0 
6 6 

= uj- 1 ^ H + T (9i, P )f{^ P ■ u-^ftE^ (2.8) 

i=l p=0 
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with Hi,E aip G Aq. We shall now specify the maps u, uj, r in detail, construct the 
corresponding M-operator and show that the Lax equation holds upon the validity of the 
GVequation of motion. 

2.2.1 Reduction of the root systems 

Let us precisely see how the root systems are embedded into each other as Aq 2 13 A# 3 D 
Aa 6 • We label the 42 roots of ^-are as 



i\p 





1 


2 


3 


4 


1 




0=2 + «3 


a=4 + «5 


a= 6 




2 


—a.2 


a=i + o= 2 + «3 


0=2+03+04+0=5 


a4 + «5 + "6 


-a 5 


3 


0=3 


0=1 + Q= 2 + 0=3 + «4 + 0=5 


«2 + «3 + «4 + «5 + "6 


a=4 




4 


— 0=4 


0=3+0=4 + 0=5 


Oi+a 2 +03+04+o 5 +o 6 


0=2 + O3 + O4 


-a 3 


5 


0=5 


0=3 + 0=4 + 0=5 + «6 


0=1+0=2+03 + "4 


Q=2 




6 


-«6 


a 5 + a 6 


0=3+0:4 


Q=l + Q=2 


— a=i 



The roots 6j, p of the A^ 6 -root system. 

We did not report values of p > 4, i.e. powers of the Coxeter element, for those roots 
which can be obtained simply by a multiplication with —1 from a root of another orbit of 
the Coxeter element. For instance 0=1,4 = — 06, 1, 0=1,5 = — 0=6,2, etc. Let us now specify the 
action of the folding map uj, which acts on the simple roots of Aq 



a; 



on for i = 1,2,3 

aj-i for i = 4,5, 6. 



(2. 



Comparing the entire root system resulting in this manner with the (£ = 3) x (h = 6) = 18 



roots of B 



(i) 



i\p 





1 


2 


3 


4 


5 


1 


6=1 


a.2 + 2q=3 


0=1 + 6=2 


— Q=l 


-(a 2 + 26= 3 ) 


— (61 + 62) 


2 


—a.2 


«1 + 6=2 + 26=3 


6=1 + 2o=2 + 2o=3 


6=2 


— (o=l + 6=2 + 26=3) 


-(61 + 26 2 + 26 3 ) 


3 


63 


0=1+0=2 + 63 


Q=2 + 63 


-6=3 


— (o=i + 6= 2 + 63) 


-(6 2 + 63) 



The roots 6j, p of the A# 3 -root system. 



it is easy to see that the map ( ^9| ) reduces the ^-root system to two copies of a I?3-root 
system plus 6 additional roots. We marked the 3 positive roots in the table of Aa & , which 
are not mapped to A# 3 via uj in bold and underlined the roots which are mapped to short 
roots in Ab 3 - The unmarked roots are therefore mapped to long roots. 

Having specified the map uj : Aa 6 — > A# 3 acting on the roots, it is important to see 
how this reduction is translated to the action on the dynamical variables. For this purpose 
we construct its "inverse" map w _1 : A# 3 — > Aa 6 , defined via the inner product relation 

u(ai)-aj = a i -uj- 1 (a j ), 1 < i < 6, 1 < j < 3. (2.10) 
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It is easy to verify that this is guaranteed by the map 

6t{ i-> uj~ l (a.i) = cti + «7_i, 1 < i < 3, ( 2 -ll) 

when taking the conventions a 2 = 2, af = a\ = 2 and a| = 1. Now we may utilize this 
map to compute the reduction map when acting on the dynamical variables q,p 

q -> u~ l {q) = uj' 1 (y^^yi&i) = (yi,V2 - yi,y3 - V2,o,y 2 - 2/3, yi - 2/2,-2/1), (2.12) 

where we used the aforementioned Euclidean realization for the A$ root system in M 7 . To 
make contact with the literature, we defined a further set of variables through the relation 
yi = ~^2 k 1 Qk, such that ( |2.12| ) becomes 

<1 (Qi, Q2, Q3, 0,-q 3 ,-q 2 ,-qi). (2.13) 



This is the reduction map as employed in [32] (see also Q). 

Likewise, we reduce next the i?3-root system to the GVroot system by means of the 
map w : Ag 3 -> Aq 2 

J ocx for i = 1,3 
a % * u(ai) = i (2.14) 
I a 2 tor i = 2. 

The "inverse" a) -1 : Ag 2 — > Ag 3 is obtained similarly as before, but now demanding 

^(a^-Qj =Q i -w~ 1 (a i ), 1 < « < 3, 1 < j < 2. (2.15) 

We find 

a\ i— > cD _1 (qi) = &i + 2^3 and «2 ► <i _1 (a2) = 3d2, (2-16) 

with the additional conventions Sf = 2 and a| = 6. The reduction map, when acting on 
the dynamical variables q,p, is now evaluated as 

q -> uj' 1 ^) = UJ' 1 (Y^. =1 Vi®i) = (m, 3m - Vi,2yx - 3y 2 ) = ?2> ft), (2-17) 

where we realized the i?3-roots in M 3 as a± = e± — e 2 , a 2 = s 2 — £3 and a.3 = £3. The 
introduction of the variables q i translates into the usual G 2 constraint q~\ + q 2 + 53 = 0, 
which corresponds to considering the three particle system in the center of mass frame. In 
fact, it will be convenient to introduce yet another set of variables, namely q^ = q 2 — (73, 
q 2 = <13 — qx an d % = q~i — q 2 to make proper contact with ( |2.1| ). 

Let us now see how to utilize these maps in order to reduce the corresponding poten- 
tials. 

2.2.2 Reduction of the potentials 

Our starting point is the ^-potential term in the form 

6 



(j 6 

VaM = -Y J ^2slV(a, p -q). (2.18) 



i=l p=0 
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Below we confirm the known fact that the Lax equation dictates that all coupling constants 
have to be taken to be the same, i.e. g^ p = g. Then Va 6 is mapped into a ^-potential of 
the form 



6 6 



Vb 3 {q) = ^J2J2T(g vp ) 2 V(uja^ 



=1 p=0 



El 

2 



E 

p=0 



2 

E 

i=l 



V(a itP ■ q) + 2V(a 3tP ■ q) 



(2.19) 



by means of the reduction map uj as specified above ( |2.9| ) and the map r acting on the 
coupling constants as 



9%v ^ t{q, 





9 

9 
K.V2 



for 5l,5 , 53,5 , 55,5 , 52,2 , 54,2 , 56,2 

for 51,1 , 91,2 , 92,1 , 52,3 , 53,0 , 53,3 , 54,0 , 54,4 , 55,4 , 55,6 , 56,5 , 56,( 

otherwise. 



(2.20) 



The map r serves here to eliminate the aforementioned additional six roots of A§ which 
have no counterpart in B% and at the same time it establishes a relationship between the 
coupling constants depending on whether the potential involves roots which are mapped 
to long or short roots. This relation is dictated by the Lax pair construction, and coincides 
with the one found by Olshanetsky and Perelomov in UJ, or [| p. 181. Note that taking 
merely the invariance of the Coxeter transformation as a guiding principle one could choose 
the coupling constants in front of the term involving long or short roots to be independent, 
see e.g. discussion in p7|]. However, integr ability demands the dependence of the coupling 
constants to be as stated in ( |2.20| ), such that one has only one coupling constant at ones 
disposal for the I?£-theories. 

Next we map the £?3-potential to the G2-potential with the help of Cj and find 



Vg 2 = ^t 



~2 5 

9 \- 



2 ^ 

p=0 



q) + 2V(ua 3>p ■ q) 



~2 

^Y,i^(^p-Q) + V(^ P -q)]. (2.21) 



p=0 



In this last reduction step we did not need to specify any additional map acting on the 
coupling constants as the embedding is now on-to-one. For consistency we re-named, 
however, g to g. 

2.2.3 Constraints from the Lax operator 

Having convinced ourselves that the potentials can be reduce properly, we still have to 



establish that the Lax operator exists and is indeed of the form (2.7). We commence by 



explicitly solving the constraint Ql.7| ) for Aq. In principle there are now 42 x 42 possible 
structure constants £%,p,j,q, with 420 of them non- vanishing. We only report here our 
conventions for the signs of the 35 essentials and leave it to the reader to obtain the 



remainin 


5 ones by 


means 


of the equations (|23 


\) and (IT 


1) 






£ l,0,l,l = 


£l,0,2,2 = 


£1,0,3,2 


= £1,1,1,2 = £1,1,2,3 = 


£1,1,3,3 = 


£1,2,1,3 


= £l,2,2,4 


= £l,3,l,4 = 


£1,3,2,5 = 


£1,3,3,5 = 


£1,4,1,5 


= £1,4,2,6 = £1,5,1,6 = 


£1,5,2,0 = 


£1,5,3,0 


= £1,6,1,0 


= £l,6,2,l = 


£1,6,3,1 = 


£2,0,1,1 = 


£2,0,2,2 


= £2,1,1,2 = £2,1,2,3 = 


£2,2,1,3 = 


£2,2,2,4 


= £2,3,1,4 


= £2,3,2,5 = 


£2,4,1,5 = 


£2,4,2,6 = 


£2,5,1,6 


= £2,5,2,0 = £2,6,1,0 = 


£2,6,2,1 = 


£3,0,1,2 


= £3,2,1,4 


= 1. 



(2.22) 
(2.23) 
(2.24) 
(2.25) 
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In fact, we verified that these choices coincide with the constants obtained directly from 
fll.3j ) when using the vector representation of A§. With fl2,22| )-( p,25| ) and the above men- 
tioned realization for the simple roots, the constraint ( |1.7| ) may be solved by 

7 

m ? a =dJ2 Va * ( 2 - 26 ) 

k=l 

which is known for some time @. Next we may solve (|L 
find 



u- 1 



for the reduced systems and 
(«)) (2-27) 



with Tfej = 1 except for T4j = 2, r 



fe=i 
k^i 

(8— i)» = 



and 



7?) 



G 2 



7 



fc=l 



(&) 



■'(ft))- 



(2.28) 



Having presented explicit solutions to the equation (1.7), we have established the existence 
of the operators L and M. In particular (2.7) and the corresponding equation for M satisfy 
the Lax equation up to the validity of the G 2 equations of motion. 

2.2.4 Conserved Charges 

It is instructive to consider an explicit matrix representation for the L-operator. Using the 
standard vector representation of Aq it follows directly form Q2.7| ) 



L 



71 



7(931) 



( ^32 /(512) 

/(?2l) ^P31 / (923,l) 

/(913) /(9l,23) ^12 

\/2/(923) V2f(m) V2f(q 21 ) 

/(92,i3) /(923) 

7(912,3) 7(^23) 

V 7(912,3) 7(92,13) 



V2f(q 32 
V2f(q 31 

V2f(q 12 


V2f(q 21 

V2f(qi 3 ) 7(91,23 



7(913,2) 7(93,12) 

7fe) 7(93,12) 

_ 7(g 3 2) 7(913,2) 

^27(912) ^7(931) V2f(q 32 



^■P21 7(923,l) 7(93l) 
^Pl3 7(912) 
V^7(923 ) 7(913) 7(92l) ^P2 3 J 



(2.29) 



where we abbreviated := Pi~Pj, qij ■= cji-qj, q~ij,k ■= qi+q~j-2qk and q kji j := 2q k -qi~qj. 
By simple matrix multiplication we compute from this the integrals of motion of the form 
I k = tr(L k )/k 

1 



h 
h 
I12 



0. 



hit 



h = n, 
5 



2-*6 



96 



-I. 4 



2 • 



h 
h 



0. 

= 0, 



:4h + 



1 



j2 



19 



i 2 > 



/ 4 = 

/10 = 

hs 



T 2 
= 0,... 



h = 0, 



1 

20 



2 • 



/ 6 / 0, h 

hi = 0, 



0, 



(2.30) 
(2.31) 
(2.32) 



12 z u 2 576 

Thus we find non-vanishing and algebraically independent charges h only for k being a 
degree of G 2 , that is 2 and 6, see e.g. [Q. Computing the L-operator of the i?3-Calogero- 
Moser model we verify the same property. 
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3. Conclusions 

We have constructed a simple expression for the L-operator of the GVCalogero-Moser 
model. We established that the constraint (1.7) may indeed be solved and therefore that 



the L and M-operator do exist. The operators are expanded in terms of H,E a £ Aq. 
The Lax equation constructed from these operators, with coefficients subject to the stated 
reduction maps, holds up to the validity of the G2 equations of motion. To find a solution 



to (1.7) and thus guaranteeing the integrability of the model, we found that we are only 
permitted to have one coupling constant in the GVtheory, instead of two, what might 
be expected from demanding invariance under the Coxeter group. In our approach this 



feature is inherited from the L?3-theory. Such a behaviour was also observed in [26]. It 
would be interesting to investigate if this limitation can be overcome by other techniques 
or to establish that this is really an intrinsic feature of the model. Furthermore, from 
the explicit computations of numerous integrals of motion, we found that they are only 
algebraically independent and non-vanishing if their degree is equal to the degree of the 
corresponding Lie algebra. 

In our discussion we did not appeal to the explicit form of the potential and only 
require the relation ( |1.6|) to be satisfied. This means that the models covered here are of 
the general form V{x) ~ 1/ sn 2 (x) including a dependence of spectral parameter fi. It would 
be interesting to investigate the properties of the spectral curves R(k, /x) = det [kl — L(/x)] 



in the spirit of [23, 24, 25, Pq| , resulting from the L-operator presented here. 

Acknowledgments. A.F. is grateful to the Departamento de Matematica, F.C.T. Uni- 
versidade do Algarve, for extremely kind hospitality and to Christian Korff for many useful 
discussions. This work was supported by the FCT project POCI/MAT/58452/2004. 
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